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NOMENCLATURE 


A 

B 

d 

M^(RN) 

q 

s 

t 

u 

V 

w 

X 

a 


Y 

6 


e 


5 


T 


C 

n(v) 

e 


plant matrix 
control matrix 
localizer beam error, m 
engine RPM, X 

pitch rate increment, rad/s 
distance from the sliding surface, s 
time, s 

surge velocity in body frame, m/s 
inertial velocity increment, m/s 
heave velocity in body frame, m/s 
state vector 

angle of attack increment, rad 
flight path angle increment, rad 

elevator angle (longitudinal stick) increment, deg (in) 
throttle increment, deg 
damping ratio 

nozzle angle increment, rad 
pitch angle increment, rad 


1. Introduction 


This report sumaarlzes the current status of our research on the 
application of Variable Structure System (VSS) theory tc design aircraft 
flight control systems. Two aircraft types are currently being Investi- 
gated: the Augmentor Wing Jet STOL Research Aircraft (AWJSRA), and 

AV-8A Harrier. The AWJSRA design considers automatic control of longi- 
tudinal dynamics during the landing phase. The main task for the AWJSRA 
Is to design an automatic landing system that captures and tracks a 
localizer beam. The control task for the AV-8A Is to track velocity 
commands In a hovering flight configuration. Much of the effort since 
our last report [1] has been devoted to developing computer programs 
that are needed to carry out VSS design In a multivariable frame work, 
and In becoming familiar with the dynamics and control problems associated 
with the aircraft types under Investigation. Numerous VSS design 
schemes were explored, particularly for the AWJSRA. The approaches pre- 
sented here are the ones that appear to be the best silted for these 
aircraft types. Examples are given of the numerical results currently 
being generated. A brief summary of VSS theory was presented In [1]. 
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2. Glide Slope Control For The 
Augnentor Wing Jet STOL Research Aircraft (AWJSRA) 


2.1 The System 

The AWJSRA Is a research aircraft modified from the De Havllland 
C-8-A turboprop by modifying the wing to Include an augmentor flap system, 
boundary layer control and other lift augmentation systems, and by 
replacing the turboprop engine by a split flow jet engine. The system 
has been described In [2, 3, 4]. The purpose of the present work was 
to design a precise glide slope control system Invariant to changes In 
some of the aircraft parameters. 


2.2 Control System Design 


The control system was designed to use the existing controls: 
elevator angle and the engine thrust - an Independent system controlled 
by the throttle. The Jet nozzle was set to a nominal value of 90°. 

The equilibrium trajectory was chosen to be a 7.5" glide at 30.9 m/sec 
(60 knots) starting at an altitude of 396.5 m (1300 ft). Other system 
parameters are given in [3,4]. The design was based on the linearized 
model In [4] modified to a wind axis coordinate system: 
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At ch« •qiiiUbrluffl trajactory, the nozzla la pazpandlcular to the aircraft's 
loaglrudlnal azla. Thus nozzle angle variation la a rather poor velocity 
control, as can be seen from the aaiall value of the control derivative. 

It was, therefore, decided to keep the nozzle angle constant. 

The variable structure control requires controls capable of almost 
Instantaneous changes. The only fast control available in the AHJSRA 
Is the elevator angle. The engine thnist control has a time constant of 
about 1 sec. The control system consis ts of two loops: the first loop 

is an internal loop with state variables [9, a-, v, qj controlled by 
the elevator angle 5^ and designed as a variable structure control system 
[5, 6]. This Is mainly an attitude control system. Speed control Is 
achieved throu^ changes In angle of attack. The control system parameters 
were determined by placement of the eigenvalues at the desired position 
when the system Is In sliding mode along the surface 

3 - (100 9) + C 2 (100 a) + C3 V + (100 q) - 0 

The eigenvalues were placed so that the resulting inner loop 
have natural frequency ot 1.5 rad/sec, damping ratio c * 0.7 
eigenvalue of 0.1 sec*^. The resulting sliding surface Is 
3 3.82(100 9) - 2.22 (100 a) - 0.934 v + (100 q) - 0 

and the resulting closed loop system, in sliding mode Is: 

6N, (2 . 4) 

h 

with q determined from (2.3); 

lOOq - -3.82 (100 9) + 2.22 (100 o) + 0.934 v (2.5) 

The velocity control through attitude can be seen by considering the 
changes In 9, a, and q required to counter a change In v, so that s->0. 

As can be seen from (2.3), a positive v can be obtained by a decrease 
In a, and vice versa. 

The use of a variable structure control system In sliding rode for 
the attitude control makes the aircraft control invariant to changes In 
the coefficients of the state matrix governing the pitch rate, q. 
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system will 
and a real 

(2.3) 


Th€t« co«fflclMits d«p«ad on the position of chs airplans e.g*. Thust 
ths AWJSRA control will bs Invariant to changas In the c.g position. 

Such a feature oay be laportant for future applications of the variable 
structure control to airplanes with large c.g, variation such as transport 
and nllltary flying vehicles. 

The outer loop control system consists of the "beam error" control 
effected through thrust variation. As the thrust direction in steady 
state is practically perpendicular to the flight path, changes in thrust 
cause changes la the vertical acceleration and, thus induce changes In 
angle of attack. As a result, a change in flight path angle occurs. 

Since beam error is proportional to y • this error is ellndnated after 
a transient motion. Changes In a also cause changes in aircraft attitude 
and velocity. These are controlled by the inner loop, which is In 
sliding mode along the surface s. The engine rpm, which controls the 
thrust l5 in turn controlled by the throttle and is unaffected by other 
state variables. 


The outer loop was designed under the assumption that the inner loop 
is already In the sliding mode. Thus, the variable q was eliminated 

using (2,5). The new state vector is thus 

- [d, 100 8, 100 a, v, 6N^] (2.6) 

• • • 

and the eq\iations of motion for 100 0, 100 a, v are those of (2.4). 


The outer control loop was designed by minimizing the quadratic 
performance index 


J - 1/2 ; [x^ Q X, + 6^1 dt 
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and is the time at which the sliding mode begins. 
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Th« rasultlag throttle control is: 


‘th - - <=i *1 (2-“) 

with 

- (5.53» 5.27, -5.91, -A. 8] (2.9) 

2.3 Sliding Motion Results 

The system dynemlcs described In (2.2) were slmuleted In the sliding 
mode, using s second order Runge-Kutts asthod. The system ves required 
to decrease a 10m. Initial beam error. The results are shown in Figures 2.1 
to 2.6. beam error decreased to SZ of its initial value in about 11 

seconds (see Fig. 2.1). There was practically no overshoot. The motion 
towards the equilibrium glide path was accompanied by a slight nose down 
pitching (.Fig. 2.2) and a very small increase in velocity ( Fig. 2.4). The 
main effect was a considerable increase in thrust as can be seen from the 
increase in r.p.m. (Fig. 2.5). This is the main path control and according 
to [2,3] is preferred by human pilots. The nose down tilt required to hold 
the speed approximately constant is also described in [3,4]. Such a control 
technique, which couples speed, attitude, and path control may be a 
heavy burden on the human pilot and, therefore, degrade his rating 
(opinion) of the system. However, the automatic control system is 
fully capable of both path, attitude, and speed control regardless of 
their coupling. 



FIGURE 2.1 The Variation of the Beam Error with Time. 

(The attitude control system in sliding mode) 



FIGURE 2.2. Variation In 6 Due to Beaa Error Control 
(The attitude control ayptem lo alldlng aoda) 




HCtURE 2.3 Variation in Angle of Attack Due to Bean Error Control 
(The attitude control sygten in sliding caode) 




RE 2.4 Variation In Velocity 
attitude control system In si 



FIGURE 2.5 Engine RPM Variation Required For Beam Error Control. 

(The attitude control system in sliding mode) 








2.4 Dftslgo of Controls Required for Reechlng the Sliding Surface 

The design procedure for reaching the sliding surface is described 
in [1, 5]. The requirement is that the norm of s should always decrease* 
l.e. * 

ss < 0 (2.10) 

The procedure described in [1] was generally used, however, since the 

rpm factor influences the motion of Che inner loop variables, an 

additional component had to be added to the control, as shown in Chapter VIII 

of [5]. From (2.1) and (2.3). we obtain 
6 

ss ■ s[ 2 a. X. - 1.2 5_] (2.11) 

1-1 ^ ^ E 

where 

a, • 0 a, « 2.14 

1 4 

a^ - .0147 aj - .173 (2.12) 

a^ • 1*03 a^ - .8 


To satisfy (2.10), the following control structure is appropriate: 
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> aj./1.2 , < a^/1.2 


(2.14) 


The following selection was made: 
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06 - 0 


(2.15) 
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The aircraft and control system dynamics both off and on the sliding 
surface were simulated and the resiilts are shown in Figures 2.7 •* 2.13. 

The response is initially slower due to a throttle command limit that 
was imposed such that 6^ < 10.7®, which corresponds to • 98.5 (normal 
take-off power setting). The general motion Is similar to that of Ideal 
sliding, with the exception that 5^ is no longer continuous (Fig. 2.13). 

The beam error (Fig. 2.7) settles In 14 seconds and Is accompanied by 
small changes In a (Fig. 2.8), nose down pitching (Fig. 2.9) and minor 
velocity variations (Fig. 2.10). The engines BPM response is given in 
Fig. 2.11, and exhibits the effect of limiting 6^. The value of s during 
the run is shown in Fig. 2.12. Note that sliding occurs almost immediately 
and is maintained throughout the maneuver. 
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3. VSS Design For The AV-8A 


This portion of our research considers the VSS design of a velocity 
coDmand control system for the AV*8A in hovering flight. Both longitudinal 
and lateral dynamics will be considered, however, this report will onlv 
addresb control of the longitudinal velocity components (surge and heave). 

Reference |6] gives the linearized model for the Harrier dynamics 
for airspeeds between 0 and 120 knots. We have elected to use the values 
for 30 knots, which result in the following model for the system dynamics 
in the body frame: 


X ■ Ax + Bu 
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In the design nozzle angle is held fixed (*^■0), so that the only n^eans 
of achieving a is by pitching the aircraft. Vertical velocity is 
controlled bv RN^. In the design of systems with a command input* it is 
customary to redefine the state and control perturbations about a commanded 
equilibrium state and control obtained by setting u ■ w - and 
solving for the remaini.ig states and controls by equating (3.1) to zero. 
This detail is omitted here but the definitions are implied. 
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3.1 VSS Design of Che ACClCude Loop 

The VSS design for sccitude control is based on Che controller 
structure shown in Figure 3.1. The sliding surface is defined by: 

8 - 56 + q , 59 - 6-9^ (3.6) 

where 

9 ■ k, 5u • 5u ■ n-u (3.7) 

Cl c 

and is regarded as a constant or slowly varying input. In sliding 
mode (s*o) we have from (3.4) and (3.6) 

9 - -6 + 9^ , - 1/Cj^ (3.8) 

which Is stable for any > 0. The transient response is dictated by 

and is invariant with respect to remaining state variables. The 
design of and is based on Fig. 3.2. The closed loop poles were 
chosen from Fig. 3.3 taken from [7]. Selecting ■ 2 rad/s and a dampine 
parameter of 3 sec ^ (which corresponds to C - .75), the resulting values 
for and k^ are: 

- 3.0 s'^ ; kj^ - .136 s/o (3.9) 

The heaving motion is controlled using a conventional proportional control 
lav 


RN - .1378 w - .8 RN (3.10) 

c 

To guarantee reaching and existance of the sliding mode, it is 
sufficient that 


ss < 0 

Differentiating (3.6) and assuming u ■ 0, we obtain 

c 

5 

ss - s( I a X. + .2 5 ] 
i-1 ^ ^ * 


(3.11) 


(3.12) 


where 
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VTOL DYNAMIC STABILITY 



Figure 3.3 Damping and Frequency Banos for Opt 
Handling Quality In VTOL Hovering Flight. 



- .0056 + C^kj^ (0.35) 
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a^ ■ Cj^kj^ (9.8) 0 


To satisfy (3.11), the following control structure Is used 
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6 - - I ’I'l *4 - s 5 ’l'4 - { 

e 1 1 s 


®1 * s > o 


8^ , s < o 
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(3.14) 


(3.15) 


Allowing for possible variations In the parameters In (3.4) with the 
exception of A(l,3), the following selections were made 
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The last term in (3.14) basically controls the time required to reach 
the sliding surface, and k^ was chosen as 4.0 in-s. 


3.2 Comparison to a Conventional Design 


A conventional design of the attitude loop for the control structure 

of Figure 3.2 follows the same lines except that the transfer function 

for 6 -► 9 is second order 
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Placing the closed loop poles to match the response time of the variable 
structure system we obtain 

■ .0816 s/m 

k 2 • 50 In 

- 24.35 in-s^ (3.19) 

Note that the gains in (3al9) are considerably higher than the gains in 
(3.14) for the variable structure control in the vicinity of the sliding 
surface (s-o) . This should aid in avoiding control saturation and 
instability due to large command inputs. 

3.3 Numerical Results 

The numerical results of this section compare the VSS control to a 
conventional control design for response stability under the presence of 
saturating control. The magnitude of 6^ was limited to 4 inches. Two 
levels of responses are given, corresponding to initial velocity errors 
of -3 m/s and -10 m/s. Figures 3.4 to 3.7 give the VSS response for an 
initial velocity error of -3 ra/s. Note from Fig. 3.6 that sliding ini- 
tiates at 2.3 seconds, when the response is essentially complete. Figure 

3.7 shows that there is little coupling with the heave d\Tiamics. Figures 
3.3 to 3.10 give the velocity pitch attitude, and longitudinal stick 
responses for an initial velocity error of -10 m/s. Note the similarity 
of response in velocity with that of Fig. 3.4. 

Figures 3.11 - 3.14 show the velocity and longitudinal stick responses 
with proportional control for the same conditions. Note that for an 
initial velocity error of -10 m/s, the proportional control is on the 
verge of instability, exhibiting 255! overshoot and prolonged periods of 
control saturation. Figures 3.13 and 3.14 should be compared with Fig. 

3.8 and Fig. 3.10. 
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TIME (s) 

Figure 3.5 Pitch Variation with VS Control for m(o) - -3 m/s 
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th VS Control for m(o) 



Surge Velocity Variation with VS Control for y(o) - -10 a/s 











gltudlnal Stick Variation 



Figure 3.13 Surge Veloclt) 




Lcngltudlnal Stick Variation With Convention Control for 



4. Future Research 


The research for the next reporting period will examine the behavior 
of the SWJSRA and the AV-8A subject to system parameter variations and 
external disturbances. In addition, transient responses will be generated 
using notillnear models for these aircraft. In the case of the AV-8A, 
we propose to examine using nozzle angle as a control to achieve reaching 
of the sliding surface. Currently, reaching takes up most of the 
transient response, and Increasing In the controller design leads to 
unstable behavior in the presence of large conmand inputs. 

For next year, we propose to examine other aircraft types currently 
of Ifiterest to NASA Ames. In particular, a tall-sltter vehicle has 
been discussed with the technical project monitor. 
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